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A BeppoSAX observation of MKN6 
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Abstract. We have used the BeppoSAX satellite to study the broad band (0.5-100 keV) X-ray spectrum of the 
Seyfert 1.5 galaxy MKN6. The source is characterized by a power law of r=1.7to'o7 an d there is no strong 
evidence for either a reflection bump or a high energy cut-off. We have detected a narrow iron line at 6.4 keV 
(rest frame) with an equivalent width of 98lgg eV. MKN6 also exhibits strong and complex absorption. At 
least two components (N// 1 = 1.341 ' 4 x 10 cm and Nh 2 =4.18I 1 ' 3 x 10 cm ) are present and they both 
partially cover the source with covering fractions of ~90% and ~50% respectively. Comparison with a previous 
ASCA observation indicates that in both absorbing columns the N_h is variable over a 2 year timescale, while 
the covering fractions are constant over the same amount of time. The state of each absorber is cold or mildly 
photoionized. The Broad Line Region (BLR) is suggested as the possible location for this complex absorption. 
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1. Introduction 

MKN6 (IC 450) is a Seyfert 1.5 galaxy which has not 
been observed often in the X-ray regime despite its bright- 
ness. The only X-ray study of MKN6 is that of ASCA 
(Feldmeier et al. 1999) which revealed a very interest- 
ing spectrum in the 0.6-9.5 keV band, showing heavy and 
' complex intrinsic absorption in the AGN nucleus, a power 
law continuum of T ~1.6 and an apparently broad 6.4 
keV iron Ka line. On the other hand MKN6 is the best 
known example of an NGC4151 analogue: both are Seyfert 
1.5 galaxies characterized by ionization cones, strong and 
complex X-ray absorption and a flat high energy spec- 
trum. All these features are quite atypical of type l objects 
and more commonly found in Seyfert 2s. Observationally 
Seyfert 1.5s are widely recognized as a distinct class of 
Seyfert galaxies (Osterbrock & Koski 1976, Cohen 1983) 
although their nature has not yet been fully understood. 
According to Nagao, Taniguchi & Murayama (2000), it 
could be that they are seen from an intermediate viewing 
angle between Seyfert Is and Seyfert 2s so that a signifi- 
cant part of the broad line region (BLR) is obscured by a 
dusty torus. Alternatively they may be Seyfert Is where 
a large part of the BLR emission is obscured by dense, 
clumpy gas clouds. Finally, they could be Seyfert 2 galax- 
ies where part of the BLR emission can be seen through 
some optically thin regions of the dusty torus. The latter 



idea is consistent with the model suggested by Feldmeier 
et al. (1999) to explain the heavy absorption measured 
in MKN6 by ASCA i.e. a dusty torus with a low density 
atmosphere. In this scenario our line of sight skims the 
surface of the torus and thus passes through the atmo- 
sphere which is assumed to have a low column density of 
the order of 10 20 -10 21 cm -2 of neutral or low ionization 
material. However the absorption measured by ASCA in 
MKN6 (Feldmeier et al. 1999) is much stronger than that 



expected in a torus atmosphere (i.e. 
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over an order of magnitude higher than that expected 
based on observations at longer wavelengths. This incon- 
sistency was overcome by Feldmeier et al. (1999) by as- 
suming either an additional absorber located inside the 
torus atmosphere or that the atmosphere is composed of 
relatively dust-free gas. 

In order to assess the true slope and the broad band X- 
ray spectrum of MKN6, to compare its characteristics with 
those of NGC4151 and to understand the nature of these 
two peculiar Seyferts, we performed an observation over 
the 0.5-100 keV band with the BeppoSAX satellite. 



2. Spectral analysis 

2.1. The BeppoSAX Observation and Data Reduction 

MKN6 was targeted by the BeppoSAX NFI in 1999 from 
September 14th to September 17th. The effective exposure 
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Fig. 1. Data to power law plus Galactic absorption model ra- 
tio. 



times were 48 ks for the LECS, 109 ks for the MECS23 
and 52 ks for the PDS. Spectra were extracted from a re- 
gion centered on MKN6 with a radius of 4 arcmin. LECS 
and MECS background spectra were extracted from blank 
sky fields using regions of the same size in detector coor- 
dinates. The net count rate was 5.02 ±0.1 x 10~ 2 cts/s in 
the (0.5-4 keV) LECS energy band, 0.25 ± 0.002 cts/s in 
the (2-10 keV) MECS energy range, and 0.5±0.02 cts/s in 
the (20-100 keV) PDS band. Source plus background light 
curves did not show any significant variation, therefore a 
timing analysis has not been performed and all data have 
been grouped for spectral analysis. 

Standard data reduction was performed using the software 
package "SAXDAS" (see http://www.sdc.asi.it/software 
and the Cookbook for BeppoSAX NFI spectral analy- 
sis, Fiore, Guainazzi & Grandi 1998). Spectral fits were 
performed using the XSPEC 11.0.1 software package and 
public response matrices as from the 1998 November is- 
sue. The spectra were rebinned in order to have at least 20 
counts per channel, this allows the use of the \ 2 method in 
determining the best fit parameters, since the distribution 
in each channel can be considered Gaussian. Constant fac- 
tors have been introduced in the fitting models in order to 
take into account the inter-calibration systematic uncer- 
tainties between instruments (Fiore, Guainazzi & Grandi 
1998). All the quoted errors correspond to 90% confidence 
level for one interesting parameter (A% 2 = 2.71). All the 
models used in what follows contain an additional term 
to allow for the absorption of X-rays due to our galaxy, 
which in the direction of MKN6 is N# Ga! = 6.4 x 10 20 
atoms cm -2 (based on 21-cm radio observations, provided 
by XSPEC). 



2.2. Data Analysis 

The broad band (0.5-100 keV) spectrum of MKN6 was 
first modeled with an absorbed power law. This model 



provides a flat photon index (r « 1.5) and an Ng w 2 x 
10 22 cm -2 , but gives a reduced x 2 of 228 for 105 degrees 
of freedom and hence is unacceptable. In fact as shown 
in figure 1, it leaves systematic residuals throughout the 
whole spectrum consisting of: (a) an evident excess emis- 
sion at low energies; (b) a marginal line-like feature at ~ 
6.4 keV; and (c) possibly a high energy cut-off. Adding 
a gaussian line to this model (model 1 in table 1), gives 
a Ax 2 of 33 for three more parameters, so that the line 
is statistically significant at >99% confidence using the 
F-test. The line is centered at 6.39^;^ keV (rest frame), 
has an equivalent width EW=114t7o eV and is consistent 
with being narrow at the 99% confidence level as demon- 
strated in figure 2; therefore in the following we have fixed 
the value of sigma to zero. The resulting flat photon index 
(r=1.52±0.03) and the count excess at soft energies could 
indicate the presence of more intrinsic absorption than as- 
sumed in the first instance. We have therefore added to 
the neutral absorber a partial covering absorption model 
(PCFABS, model 2) which has been successfully applied 
to this object in the ASCA observation (Feldmeier et al. 
1999). The addition of this partial covering absorber im- 
proves the quality of the fit at greater than >99% signif- 
icance level (Ax 2 = 103 for two more parameters). Using 
this model, the source is completely covered by cold uni- 
form material with N# ~4-9 x 10 21 atoms cm~ 2 , and 
partially covered (~ 70%) by another with N H ~3 x 10 22 
atoms cm~ 2 . Despite the goodness of this fit some resid- 
uals below 1 keV are still present. 

We have investigated if the presence of a diffuse hot gas lo- 
cated outside the intrinsic absorption and heated by star- 
burst activity in the host galaxy, could be responsible for 
the soft excess component. The addition of a thermal com- 
ponent (Raymond- Smith model in XSPEC) to model 2 
does not improve the quality of the fit (x 2 =203 with 99 
d.o.f). Further, this model was also rejected by Feldmeier 
et al. (1999) when applied to the ASCA data. 
If we substitute in model 2 the neutral absorber with an- 
other partial covering one (model 3), the fit improves fur- 
ther (>99%, Ax 2 =8 with one more parameter) and pro- 
vides our best fit model to the data as also found with 
the ASCA data. Since the width of the line is sensitive to 
the underlying continuum, we have left the width free to 
vary in model 3 and also in this case the line is consistent 
with being narrow. This is also true for any of the models 
listed in table 1. In figure 3 the broad band (0.5-100 keV) 
spectrum of MKN6 fitted with model 3 is shown. 
An alternative model used to fit a multicolumn absorber is 
the dual-absorber. This model requires different assump- 
tions on the source + absorber geometry and thus a differ- 
ent interpretation of the data. In the dual-absorber model, 
two different columns cover the source and the relative 
normalization of the power law absorbed by these two 
columns gives the percentage of the source covered by each 
Njj in such a way that the sum of both percentages is 100. 
We have applied this model to our data (model 4) but no 
improvement in the quality of the fit has been found. 
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To check the robustness of our findings, we have also inves- 
tigated the presence of a warm medium instead of a cold 
one via the ABSORI model in XSPEC. In the case when 
the cold medium is represented by a uniform absorber the 
X 2 worsens significantly (195/101) thus this model must 
be rejected. If instead the cold medium is assumed to be 
an absorber partially covering the source, a combination 
which was successfully used to model the complex ab- 
sorption of NGC4151 (Schurch and Warwick 2002, Piro 
et al. 2002), then we obtain a statistically acceptable fit 
to our data. However we are not able to constrain Nh 

± 1 w a t m. 

and £ with this model. Reasonable values of the ioniza- 
tion parameter, deduced from previous studies of Seyfert 
galaxies (Perola et al. 2002) and compatible with the lack 
of strong warm absorption features in the spectrum, all 
provide acceptable fits: for example assuming a column 
density similar to that obtained in the case of uniform 
covering, N Hmarm ~3 x 10 22 cm~ 2 and £ w 10 erg cm s _1 
gives x 2 =109/100. Therefore we cannot exclude that the 
absorbing gas is just mildly ionized on the basis of our 
data, although we are not able to constrain its properties 
significantly. 

We have also searched for the presence of a high energy 
cut-off and a reflection component which are typically ob- 
served in Seyfert 1 galaxies (e.g. Perola et al. 2002) by 
replacing the power law with the PEXRAV model (model 
5). However, we do not find strong evidence for either of 
these two components and we are only able to put upper 
limits on the parameters values: R<1.2 and E c >70 keV. 

2.3. Comparison with the ASCA results 

In this subsection we briefly compare two X-ray studies of 
MKN6: the ASCA observation performed on April 1997 
and our BeppoSAX measurement from September 1999. 
The source shows flux variability in the 2-10 keV energy 



range; in fact the ASCA absorbed (unabsorbed) flux is 
-1 (1.4) x 10" 11 erg cm" 2 s" 1 while for BeppoSAX it 
is ~ 2.5 (3.1) x 10~ n erg cm~ 2 s _1 . These flux values 
are obtained using the same model (model 3 in table 1) 
while the absorption correction is performed by forcing 
both column densities to zero. 

The BeppoSAX observation of MKN6 basically confirms 
the ASCA results i.e. the detection of heavy and complex 
X-ray absorption towards the source nucleus. In both ob- 
servations a double partial covering model provides the 
best statistical fit to the data but from the comparison 
of our results with those of ASCA a variation in the col- 
umn densities has been detected. Our analysis of ASCA 
data provides N Hl = 2.99±qH x 10 22 cm" 2 with a cov- 
ering fraction of -90% and N H . 2 = 11±|;^ x 10 22 cur 2 
with a covering fraction of —50%. These values are con- 
sistent within the errors with those reported by Feldmeier 
et al. (1999). BeppoSAX detected smaller column densi- 
ties of 1.2 and 3.6 x 10 22 cm~ 2 with covering fractions 
C fl =93±%%% and C/ 2 =52t^% respectively. Therefore, 
while the covering fractions C / x and C / 2 are not changed 
with respect to the ASCA observation, within their indi- 
vidual statistical uncertainties, a significant variation in 
the column densities of the absorbing materials is found 
over a two year timescale. This is evident in figure 4 
where the LECS+MECS BeppoSAX data are fitted with 
model 3 but imposing the best fit parameters obtained 
from ASCA: clearly the two data sets cannot be recon- 
ciled without changes in the spectral shape parameters. 
To further prove absorption variations, we show in figure 
5 the confidence contours of both column densities in the 
BeppoSAX (solid lines) and ASCA (dashed line) observa- 
tions. Since both covering fractions are very similar and 
relatively well constrained, we have fixed them to their 
best fit values. A check for variations in the absorption 
over a shorter time-scale has also been performed. Since 
the BeppoSAX observation does not show any significant 
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Fig. 4. LECS+MECS BeppoSAX data fitted with model 3 but 
imposing the best fit parameters obtained by ASCA for both 
partial covering absorbers. 
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Fig. 5. Comparison of 1 the column densities found by ASCA 
(dashed lines) and BeppoSAX (solid lines). 



variation in the continuum light curve, the whole obser- 
vation (about 100 ks MECS exposure) has been divided 
in two roughly equal parts and the same spectral analysis 
of the entire observation has been repeated on both seg- 
ments of the observation. We find no significant variation 
in any of the spectral parameters, and in particular in the 
amount of absorption. We can therefore conclude that no 
short timescale (one day) variations are present. 
Regarding the iron line, we note that the broadening seen 
by ASCA (<t ^230 eV) is still consistent with our mea- 
surement of a narrow feature (see figure 2). 
Finally, we measure a 

slightly harder spectrum (r=1.7to^o?) than 1.6 as fixed 
in the ASCA data. 



3. Discussion and Conclusions 

The MKN6 spectrum observed by BeppoSAX resembles 
that of NGC4151 in the absorption properties although 
the broad band spectral characteristics (the power law, 
reflection component and high energy cut-off) are similar 
to more classical Seyfert 1 objects (Perola et al. 2002). The 
absorption is complex and variable. At least two compo- 
nents, each with its own column density and covering frac- 
tion, are present. In both columns the N# is variable over 
a two year timescales while the covering is constant over 
the same time. The state of each absorber is cold or mildly 
photoionized, in fact a combination of cold and warm gas 
as seen in NGC4151 is not ruled out by the data. Recently 
Maiolino et al. (2001) found that the E(B-V)/N ffx ra- 
tio estimated from optical/infrared broad band lines and 
X-ray data is lower than Galactic in a sample of AGNs, 
including MKN6; this issue was discussed in detail also 
by Feldmeier at al (1999). Our results of a variable N# x , 
clearly suggest caution in adopting such a ratio as an in- 
dicator of anomalous absorption properties. Furthermore, 
the optical band lines in MKN6 are known to be variable 
in time (Sergeev et al. 1999); therefore an estimate of the 
E(B-V)/N# ratio strongly requires simultaneous optical 
and X-ray measurements. 

The variation observed in N# can be explained either by 
a change in the ionization state of each absorber due to 
a change in the incoming radiation and/or to variation in 
the amount of the absorbing gas along the line of sight. 
A joint fit of ASCA/GIS and BeppoSAX/MECS data are 
compatible with a change in the £ parameter although the 
limitation of the ABSORI model in XSPEC does not al- 
low us to firmly prove a variation in the ionization state of 
the absorption. The other possibility is that the amount of 
gas along the line of sight varied between the ASCA and 
BeppoSAX observation. In this case the absorber must 
be clumpy and the variation timescales are related to the 
typical crossing time of an absorbing cloud along the line 
of sight. Following the reasoning of Risaliti et al. (2002) 
the location of the absorber is related to the black hole 
mass, the column density of the absorber, the cloud den- 
sity and the absorption variation: the BeppoSAX data in- 
dicate that the variability timescales are greater than 1 
day and smaller than 2.4 years. For reasonable values of 
the other relevant parameters we find that the molecular 
torus and the BLR are both viable locations for the ab- 
sorption (Risaliti et al. 2002). In the first case we go back 
to the original suggestion of Feldmeier et al. (1999) who 
propose the atmosphere above the torus as the X-ray ab- 
sorber. However the column densities measured both by 
ASCA and BeppoSAX are substantially larger than ex- 
pected in a torus atmosphere (<~ 10 20 - 10 21 cm~ 2 , Wilson 
1996). Considering these difficulties, it is reasonable to 
consider the BLR as an alternative scenario. The BLR lo- 
cation would be consistent with our variability constraints 
as well as with the value of column densities and covering 
fractions observed. Also the requirement for a dust free 
gas, indicated by the greater absorptions seen in X-rays 
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N Hl x 10 22 




Nj/ 2 x 10 22 
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E Ltne (keV) 
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(3) 
(4) 
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OQ+0 02 
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1 66+ 07 

i - uu -0.06 

i 7n +0 08 

1 66+ 09 

i - uu -0.05 

1 63+ ' 08 

J - uo — 0.06 


6.39^? 

6.40+ni 
6 42 +0 - 15 

u -^-0.14 

6 42+ ' 14 
6 42+ ' 04 

u -* 4 — 0.06 


11 -70 

97t« 
137+H 

qo+34 


<1.2 


> 70 


195/103 
92/101 
84/100 
92/101 

83/98 



than at longer wavelengths, is compatible with the BLR. 
Gas lying here is quite likely to be dust free since this re- 
gion is within the sublimation radius of the central engine. 
Furthermore, the partial covering model itself requires 
that the location of the absorber is close to, and com- 
patible in size, with the emitting region thus making the 
BLR a likely site (Reichert ct al. 1986). However in this 
scenario, given the complexity of the absorber, one must 
assume that the BLR itself is either obscured by dense 
clumpy gas clouds or is itself complex. This last sugges- 
tion is in line with reverberation studies which show that 
the BLR is slightly stratified into multiple emitting zones 
with strong gradients both in the ionization parameter 
and/or density (Baldwin et al. 2003). It is therefore likely 
that changes in the central engine produce variations in 
one or more of these stratifications. It remains to be un- 
derstood why such heavy and complex X-ray absorption is 
only present in a few objects and not in all Scyfert Is. The 
fact that both NGC4151 and MKN6 are Seyfert 1.5 may 
bear some relation to the peculiarity in their absorption 
properties. Based on the results of our analysis it is likely 
that this enigmatic class is made of Seyfert 1 's with a BLR 
cither obscured by extra clumpy gas clouds or peculiar in 
its geometry. 
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